SISSA 31/2004/FM 



Form factors in the SS model 



and its RSOS restrictions 



Benedicte Ponsot 1 

International School for Advanced Studies (SISSA), 
Via Beirut 2-4, 34014 Trieste, Italy 



Abstract 

New integral representations for form factors in the two parametric SS model are pro- 
posed. Some form factors in the parafermionic sine-Gordon model and in an integrable 
perturbation of SU (2) coset conformal field theories are straightforwardly obtained by dif- 
ferent quantum group restrictions. Numerical checks on the value of the central charge are 
performed. 

PACS: ll.10.-z ll.10.Kk 11.55.Ds 

1 Introduction 

Form factors with n particles of an operator are the matrix elements of this operator between the 
vacuum and a n particles state. They are non pertubative objects that can be constructed, up to 
a normalization, as solution of "bootstrap" equations [1-3], once the exact scattering matrix is 
known. Besides being solutions of a nice mathematical problem, the form factors are useful tools 
to determine the long distance expansion of correlation functions by inserting a complete set 
of asymptotic states. It is often enough to approximate (with a good accuracy) the correlation 
functions of local operators with the contribution of the form factor with the smallest number 
of particles, due to the fast convergence of the spectral series 2 . 

In this paper, we construct form factors in a two parametric family of massive integrable 
quantum field theories known as the SS model, whose action can be found in [4], and which can 
be written in terms of three boson fields ipi, i = 1,2,3 with an exponential interaction. We will 

1 ponsot@fm.sissa.it 

2 It is usually believed that the spectral series converges for local operators, but this has not been proven so 

far. 
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restrict ourselves to the unitary regime of this theory [4], where the model has a U(l) x U(l) 
symmetry described by two conserved topological charges: 

Q± = \{Qi ± Q2), Qi = J dx 1 ^, i = l, 2 , 

where j®' 1 are the components of the U(l) current. Both classically and in the quantum case, 
the charges (eigenvalues of the charges) satisfy the conditions: 

Qi,Qi€Z, Qi + Q 2 G2Z. 

The U(l) x U(l) symmetry can be extended up to the symmetry generated by two quantum 
affine algebras U qi (sl2) <8> f/g 2 (s£ 2 ). The name of the SS model is unfortunately justified by 
the expression of its two-particles S-matrix -first considered in [5]- in terms of two sine-Gordon 
5 s ' G -matrices with different coupling constants (3\ , /3 2 : 

S(9 12 ) = S^(9 12 ) ® S^{9 12 ), 

where pi = g7r _^a , Pi = 87r _^a [% and #12 = Q\—Q 2 is the rapidity difference. The unitary regime 
is characterized by the following conditions on the parameters: p\,P2 > and p\ + P2 > 2. We 
will need to introduce the negative parameter p% such that Pi +P2 +P3 = 2, and the parameters 
aj, i = 1,2, 3, defined by Pi = 2a? [4]. 

The SS model includes, as particular cases, the N=2 supersymmetric sine-Gordon model, whose 
S'-matrix has the form [7]: 

as well as the 0(4) non linear sigma model (Principal chiral field model [8,9]), which possesses 
577(2) x 57/(2) symmetry; its 5-matrix reads [6] 

5(^i2) = -5^(^2)0 5^(^2), 

which is nothing but the tensor product of two 5-matrices of the 577(2) invariant Thirring 
model. It also contains the anisotropic chiral field [10], which has the U(l) x SU(2) symmetry; 
its 5-matrix is 

S(9 12 ) = -S^ G (9 12 )^S^ G (9 12 ). 

The 55 model also includes other known integrable QFTs, like the 0(3) non linear sigma model, 
the sausage model [11] and the cosine-cosine model [12]. Let us note in passing that the ground 
state energy of the model in finite volume was determined in [4] by TBA method, and the finite 
size effects for the 55 model were recently studied in [13]. 

In order to construct the form factors in this theory, we give in the first section a summary 
of the method used in [14, 15] to construct form factors in the sine-Gordon model, and recall 
how to perform RSOS restriction [16] directly on sine-Gordon form factors [17]. In the second 
section, we propose a generic formula inspired from the SG model for form factors containing 
an even number of particles in the 55 model; in particular we write down the form factors 
of the trace of the energy momentum tensor, first obtained in [5], and of some of the expo- 
nential fields (this problem was also recently considered in [18], with different methods). Two 
particles form factors of these operators can be expressed in terms of two particles form fac- 
tors in the SG model. Then we do two different kinds of RSOS restrictions, namely [19] to 
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the parafermionic sine-Gordon model (this model includes, amongst others, the N=l and the 
restricted N=2 supersymmetric sine-Gordon models), and to the integrable perturbed coset 
CFT [20] sit(2) pi _ 2 ® •su(2)p 2 _ 2 /sw(2)p 1+ p 2 _4. A common limit of these QFT's is the Polyakov- 
Wiegmann model [8]. We propose expressions for form factors of the trace operator (they were 
first considered in [5]). In particular, we check our expressions for the two particles form factor 
of the trace of the energy momentum tensor by making a numerical estimation of the central 
charge and compare it with the exact result. 



2 Form factors in the sine-Gordon model and RSOS restriction 

In this section we recapitulate known results on form factors in the SG model in the repulsive 
regime. They will be useful in the next section for the construction of the form factors in the 
SS model. 

The Sine-Gordon model alias the massive Thirring model is defined by the Lagrangians: 



respectively. The fermi field ip correspond to the soliton and antisoliton and the bose field ip to 
the lowest breather which is the lowest soliton antisoliton bound state. The relation between 
the coupling constants was found in [21] within the framework of perturbation theory: 

f3 2 ix 

P = 



8vr-/3 2 TT + 2g' 

The two soliton sine-Gordon S'-matrix contains the following scattering amplitudes: the two- 
soliton amplitude a p (6), the forward and backward soliton anti-soliton amplitudes b p (9) and 

c P (ey. 

smh9/p sinhivr/p 
sinh(«7r — 9)/p sinh(z7r — 9)/p 



... f°° dt sinh ^(1 — p)t sinh 

a p (6) = exp / — 

Jo * cosh \ sinh ^pt 



dt sinh |(1 -p)t sinh ^ 
| sinh If] 

This S'-matrix satisfies the Yang-Baxter equation as well as the unitarity condition: 



s^ G (e)s s p G (-e) = i, 

which can be rewritten for the amplitudes as: 

a p {9)a p {-9) = 1, b p {9)b p {-d) + c p {9)c p {-6) = 1. 
The crossing symmetry condition reads for the amplitudes: 

a p (m - 9) = b p (6), c p (m - 9) = c p (9). 
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The form factors f(0i, • • • , 02n) of a local operator 3 in the SG model are covector valued functions 
that satisfy a system of equations [3], which consist of a Riemann-Hilbert problem: 

f(Ql, ■ ■ ■ • • • 5 ^2n)Sp G (9i — Oi+l) = f(0l, ■ ■ ■ , Qi+l,0i, ■ ■ ■ , 02n), 

f(01, ■ ■ ■ , 02n-l,02n + 2m) = f(02n, #1, " ' ' > #2n-l), 

and a residue equation at 0\ = 6*2 n + in '■ 

(2n- 1 \ 
1 - IJ 5"p G (^ - 02n) J eo, e = Si ® S 2 „ + Si (g> S 2 n, 

where s (or +) corresponds to the solitonic state (highest weight state), and s (or -) to the 
antisolitonic state. Form factors containing an arbitrary number of particles for the energy 
momentum tensor, the topological current and the semi-local operator e ±l 2^sG were first con- 
structed in [3] by Smirnov; form factors of non local exponential fields e mvsG were constructed 
in [22] by Lukyanov, using free field representation techniques that provide integral representa- 
tions different from those of [3]. Below, we make the choice to present the posterior construction 
presented in [14, 15] by Babujian and Karowski et. al. 

We first introduce the minimal form factor f p {0\2) of the SG model: it satisfies the relation 

f P (0) = -f P (-0)a p (0) = f p (2i7r - 0), 

and reads explicitly 

ro ° dt sinhi(l-p)/ 1- cosh £(!-£) 



f p (0) = -isinh-/™"(#) = -*sinh-exp J 



t sinh \pt cosh \t 2 sinh t 
Its asymptotic behaviour when — ► ±oo is given by f P (0) ~ C p e ±li{ p +1){e iw \ with the constant 



C =^expl r dt ( sinh Kl-^ l-p\ (2) 

p 2 L 2J t y sinh \pt cosh \t sinh t pt J ' W 

It is proposed in [14] that form factors in SG can be generically written 4 : 

f(9 1 ,02,...,9n) =N n Y[f p (O ij ) f dm... f du m hp(0,u)p n (0,u)*P(0,u), (3) 
j<j JCe JCg 

where we introduced the scalar function (completely determined by the 5-matrix) 

2n m 

h p (0, U) = \ '{ jj (f>p(0i ~ Uj) l\ T p (u r - U s ), 
i=l j=l \<r<s<m 



with 



M u ) = f 7 w ( -— "T' T p( U ) = 



fp{u)f P (u + in) 4>p(u)4> p (-u) 

3 We shall consider in the following form factors with an even number of particles only. 
4 This representation holds whether operators are local or not, topologically neutral or not. 
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ty p (9, u) is the Bethe ansatz state covector: we first define the monodromy matrix T p as 



A{e 1 ,...,e n ,u) B(e 1 ,...,e n , u ) 

C(9 1 ,...,9 n ,u) D(9 1 ,...,9 n ,u) 



T p {e x 



...,9 n ,u) = S* G (9 1 -u)...S: 



{6n ~ u), 



the definition of the Bethe ansatz covector is given by 



i& p (6, u) = 0,1. ,. n C(6>i, . . . , 9 n ,ui) ... (7(6*1, ... , 



where fii... n is the pseudo vacuum consisting only of solitons 



.71. 



= s <g> ■ ■ • <8> s. 



The number of integration variables m is related to the topological charge Q (Q £ Z) of the 
operator considered and the number n of particles through the relation Q = n — 2m. For 
example, for n = 2 and Q = 0, 

¥ p (0i,0 2) u) = ^_( i^2,^) + ^ p + (^i,02,m) = b{9 l -u)c{9 2 -u)s 1 ®S2+c{9 l -u)a{9 2 -u)s l ®S2. 

The function p n (9,u) is the only ingredient in formula © which depends on the operator con- 
sidered. If the operator is chargeless, the form factors contain an even number of particles, and 
if in addition the operator is local, then the p- function satisfies the conditions 5 : 

1. P2n(9, u) is a polynomial in e ±u i , (j = 1, . . . , n) and P2n(9, u) = P2n(- • • , #i — 2iir, . . . , u) 

2. P2n{9\ = 92n+iTT, . . . 9 2 n,Ul ■ ■ ■ U n = 9 2n ) = P2n~2(9 2 ■ ■ ■ 9 2n -l)Ul ■ ■ ■ Un^+p 1 (9 2 . . . 2n -l), 
P2n(0l = 92n+iTT, . . . 9 2n \Ul . . . U n = 9 2n ±i^) = P2n-2(#2 • ■ -^n-l^l • ■ ■ U n _ l )+p 2 ± (9 2 . . . 02n-l), 

where p ' 2 (9 2 . . . 9 2n —\) are independent of the integration variables. 

3. p 2n (9,u) is symmetric with respect to the 0's and the it's. 

4- p 2n {9 + In A, u + In A) = A s p 2n {9, u) where s is the Lorentz spin of the operator. 

Finally, the integration contours Cg consist of several pieces for all integration variables Uj : a 
line from — oo to oo avoiding all poles such that Im0j — ir — e < Imuj < Im9i — tt, and clockwise 
oriented circles around the poles (of the 4>(9i — Uj)) at 6\ = Uj, (j = 1, . . . , m). 

Trace of the energy momentum tensor. The trace operator is a spinless and chargeless local 
operator. Its p-function is [15]: 





The residue equation gives the following relation for the normalization N® n (see [14]): 




5 We consider here only the case where the operator is of bosonic type and the particles are of fermionic type. 
If both are fermionic, there is an extra statistic factor to be taken into account, see [15]. 
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The two particles form factor can be computed explicitly: 

2ttN® cosh — 

fsaWvi) = —rr-fpWn) . . lf . 2 - > i + s 2 ), 
L£ smh ^(nr - 6 12 ) 

in agreement with the result first obtained by diagonalization of the S^-matrix in [1]. The 
normalization for two particles is chosen to be Nf = ^iM 2 Cp (M being the mass of the soliton), 
in order to have: 

/#g(0i + i7T , °i) = 2vrM 2 (si g> s 2 + Ii <8> s 2 ), 

Exponential fields. The p-function of the spinless, chargeless, non local, exponential field e ia<PSG ^ 
in the SG model is not written in [15], though it is a minor extension of the results contained 
in this paper. It reads 

2aUj 

P l G (e,u) = 4^ ul : ie 1 • (5) 



e - 







The conditions 1. and 2. are modified into: 

— 2z7ra — 

1. p| G (0,u) =e"— P § G (...,6i-2iir,...,u), 

2- Ps G (^l = #2n + ivr, . . . 6*2 n ; Ul ■■■U n = 2n ) = e P Ps G (02 ■ ■ • 02n-l] «1 • • • «n-l), 
PsG^l = ^2n + ivr, . . . 2rl ;Ul • • • U n = #2n + «0 = Psq(02 ■ ■ ■ 02n-V, Ul . . . lt„-l), 
^5g(^1 = ^2n + ivr, . . . 2n ] Ul ■ ■ ■ U n = 2n ~ «7r) = e ^ P5 G (#2 • • • 02n-ll «1 • • • 

The form factors of exponential fields were first constructed in [22] with a different representation. 
For a = with fcgZ, the form factors can be computed explicitly, and their expression can 
be found in [22]; In particular, the form factors of the semi-local operator e^a'Psaim) were fi rs t 
obtained in [3]; their expression with two particles is 6 : 

= fp{e vl )—— y- -— e ^^si®S2 + e T ^^si®S2 

q}sinh±(z7r - 12 ) V 

An important remark to be made is that there exists an alternative p e -function to eq. (j3J) for 
the trace of the energy momentum tensor. Indeed, if one remembers that in the SG model the 
trace is identified as the term in the action cos (3(psG, U P to inessential coefficients, then one 
can rewrite its p-function in a suggestive form as a sum of p- functions for exponential fields with 
a = P and a = — /?, i.e.: 

p1 G (0,u)=p 1 SG (0 1 u)+p- G (0,u). (6) 

Indeed, using the expression 

, , , , 2tt cotan^2 cosh^ / . ^-»i2 ^ ^-e 12 \ 

S smh±(>7r - (9i 2 ) V / 



We will not need the normalization factor of the exponential fields in the SG model. They can be found 
in [23]. 
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we find 

ft , . 2vriV 2 cotan^ , , cosh ^ 

fs G (0u) = ^4 ~/p ■ , i r „ M ® S "2 + S "l ® 

C-p smh ^(zvr - 6*12) 

where the new normalization constant 2V® is given this time by the formula iV® = tM p Cp tan ^ . 
Let us note that the vacuum expectation value of the trace was obtained in [24] thanks to the 
thermodynamic Bethe ansatz, and is equal to < >= — ttM 2 tan the following relation then 
holds: 

at® = - < e > X 

P 7T 

RSOS restriction [16]. The RSOS restriction describes the ^i^-perturbations of minimal mod- 
els of CFT [25] for rational values of p. When p is an integer, we deal with ^i^-perturbations 
of minimal models M p with central charge c = 1 — p ^ p+ i^ ■ I n particular, we remind that 
gRSOS = — 1 is the Ising S'-matrix. Form factors in the model M p can be directly obtained 
from those of the SG model, as explained in [17]. For the trace operator, the RSOS procedure 
consists of 'taking the half of the p-function 7 (j3J), such that its p- function reads: 

2n n 

P%sos(0^) = -T, e ~ 6 'T, eUm ' 

i=l m=l 

then we should modify the Bethe ansatz state: 

2n 
i=l 

The two particles form factor of the trace operator reads explicitly: 



vj/P = e^^^vj/P ei = ± Y^e, = (7) 



n 27T 2V® COSh ^ hl _ i7T 

C£ smhi(ivr - 12 ) 

and the normalization for two particles is chosen to be 2V® = ^iM 2 C p , in order to have: 

fRSOs(°i + «r,0i) = 2vrM 2 (et si8)s 2 + e"t s 1 ®s 2 ). 

The knowledge of the form factors of the trace of the stress energy tensor allows to estimate the 
variation of the central charge by means of the so-called "c-theorem" sum rule [26,27]: 



3 

cuv = 



/•oo 

/ drr 3 < 9(r)0(O) > . (9) 

Since in the massive case any correlation function can be represented by its spectral expansion 
< O{x)O{0) >= (10) 



7 The rationale behind this is explained in [17]. 
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V 






Q 
O 


n ^ 


n ^ 


o. 1Z 


U.OOOl 


U.tJOO^O'i 


A 


u.uyoo 


n 7 


5 


0.7972 


0.8 


10 


0.9373 


0.9454... 


20 


0.9744 


0.9857... 


100 


0.9864 


0.9994... 



Table 1: Minimal models M, 



the computation of c turns out to be a non trivial check for the form factor f^sosi^ 12 ) m ec i ua tion 
(JHJ. In Table 1 above are presented the numerical results with two particles contribution for the 
central charge in the M p model versus the theoretical result. As one may observe with the case 
p = 3.12, the parameter p can be taken continuous, as the observables depend continuously on 
it 8 . 

More generally [17], one can obtain the form factors of the primaries $x,fe using the identification 
$i,k ~ e* ( 2 -P<pso{®)^ anc [ then 'twisting' the Bethe ansatz state like in (|7|). In the two particle 
case, the form factors can be computed explicitly. 



3 Form factors in the SS model and its RSOS restrictions 

This program (for pi,p2 > 1) was initiated for this model first by Smirnov [5] for form factors of 
the energy momentum tensor and the U(l) x U{\) current, then continued recently by Fateev 
and Lashkevich [18], who proposed expressions for the form factors of the exponential fields 
using the method of [22]. 

The form-factors F{6\, • • • , 02n) of a local operator in the SS model satisfy the system of equa- 
tions [5]: 

F(9i, • • • , 9i, 0i + i, • • • , 02n)S(0i — 9i + \) = F(0i, • • • , 8i + i,8i, • • • , 02n), 

F(9 l , ■■■ , e 2 n-i,0 2 n + 2in) = -F(9 2n , 0i, ■ ■ ■ , 9 2n -i), (11) 

and the residue equation at 9\ = #2n + i?n 

(2n-l \ 
i - n s m (*< - e ^ s p* & - 02n n e °' ( 12 ) 

e = (si s 2 + s\ <g> s 2 ) <8> (si (g) s 2 + s\ <g> s 2 ). 

The minimal form factor satisfies: f S s(9) = — f ss (—9)a pi (9)ap 2 (9) = f ss (2i7r — 9), and reads [5]: 

fss(e 12 ) = —gtf Pl (0l2)fp 2 (0ll). (13) 



We refer the reader to the Fig. 5 of [28] for similar numerical tests on the central charge (the parameter p is 
taken continuous); in this article, the author constructed form factors with two kinks only, starting directly with 
the definition of the RSOS matrix. It is trivial to see that our correlation function < OO > with two particles 
coincides with the one of [28]. 
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It has no poles and no zeros in the physical strip < Im 9i 2 < vr and at most a simple zero at 

012 = 0. 

We make the following ansatz for the form factors in the SS model 9 : 
F(9i, . . . , 9 n ) = 

AA n J]/ ss (%) f du x ... I du s h Pl (9,u)^(9,u) f dv x ...f dv t h P2 (9,v)V P2 (9,v) 

j<j J Co JCe J Co JCe 

xM n (9,u,v) p n (9,u,v). (14) 

A few comments about our ansatz are in order: 

• in the equation (fT2)) . the residue at the pole located at 9\ = 9 2n +in is simple. As explained 
in [5], the introduction of the cosine term in the minimal form factor ()13|) gives a zero at 
#i = #2n + in, and consequently we should look for the residue at a second order pole of 
the integral representation proposed in equation ©. 

• the properties for the p-function of a local operator are similar to those of the SG model 
presented in the previous section, though with a small modification at 9\ = 9 2n + in: 

P2n(0l ■ ■■02n,Ul . . . U n , V\ . . . V n ) = -J>2n-2(#2 • • • #2n-l,«l • • -U n -l,Vl . . . U n _i) 

+ p 1 {9 2 . . . 9 2n -i) at u n = v n = 6> 2n , 

P2n{dl ■ ■■02n,Ul ...U n , V\ ... V n ) = Pln-2^2 • • •#2n-l,'Ul • • • U n -X,V\ . . . U„_l) 

+ £±,±(#2 • • • 02n-i) at u n = e 2n ± in, v n = 9 2n ± in. 

• clearly [5], the simple tensor product of the form factors of the SG model is not a solution 
as it spoils the residue equation. At Q\ = 9 2n + iir, each of the 2n integration contours 
get pinched at u n = 9 2n ,9 2n ± m and v n = 9 2n ,9 2n ± in (the choice for the integration 
variables u n and v n is arbitrary because of symmetry). For this reason we introduced the 
function M 2n (9,u,v) which must have the properties at 9\ = 9 2n + in: 

- M 2n (9i . . . 9 2n ,ui . . .u n ,v± . . . v n ) = M 2n - 2 (9 2 . . . 9 2n -i,u 1 . ..u n ~i,vi . . .t> n -i) 
when u n = v n = 9 2n or u n = 9 2n ± in and v n = 9 2n ± iir. 

- M 2n (9i . . . 9 2n , u\ . . . u n , vi . . . v n ) = when u n = 9 2n and v n = 9 2n ± in, or when 
v n = 2n and u n = 9 2n ± in. 

We introduce the set S = (1, . . . , 2n) as well as T C <S and T = S\T. A new result of this article 
is the following expression 10 : 



M 2n (9,u,v)= 



IlMef COS 2 ^ 



2/ 



2^i=i e 2^i=i e res, H»eT, sin -gf 



#T=n-l 



n «*■. cos^ n ier, cosfc? 

x r=l,...,n r=l,...,n 

ril<r<s<n. COS ^2^ Yll<r<s<n COS ~ 2^" 



The numbers s and t of integration variables u and v are related to the topological charges Qi and Q2 as well 
as to the number of particles n through the formula Qi = n ~ 2s and Q2 = n — 2t. 
10 The notation stands for 'number of elements'. 
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For two particles, this function is equal to M.2(8i, 9 2 , u, v ) = 1. 

The function above is not the only one to have the required properties, but we conjecture that 
this is the one we need. 



Trace of the stress energy tensor. We propose the following p-function: 
with the recursion relation for the normalization constant: 



/ frnin(n\ fminfn\\^ 
Jy 2n — ~ lJV 2n-2' 



87m 2 

The p^-function is given indifferently by eq. Q) or eq. ©. In the two particles case, the 
form factor boils down to the sum of products of two form factors of the SG model, and reads 
explicitly: 

-F (#i, #2) = — cosh^ 2 - f ss {9 l2 ) ( (Sl0S " 2 1 +g " l0S2)0(gl ^" 2 + S " l0g2) l . (16) 



It is not difficult to check that this expression satisfies equations (|11|) : this result agrees with [5]. 

The constant N 2 was chosen equal to N 2 = ^^^^r~~ ' sucn that the following relation holds: 
F®(9i + iir, 9\) = 2ttM 2 eo, where M is the mass of the fundamental particles. 
For p\ = p 2 = 1, P3 = 0, we have two sine-Gordon models at the free fermion point; the 
expression of equation ()16j) becomes as expected: 

F @ (9 1 ,9 2 ) = -2vrM 2 sinh^ 2 - e . 

A few numerical checks on the value of the central charge at the level of two particles con- 
tributions can be found in Table 2. Compared with the results of Table 1, one sees that the 
accuracy is rather poor. We have so far no compelling explanation for this phenomenon. It 
worth mentioning that in [29] , it was already noticed that the two-particle approximation to the 
correlation function does not always give accurate results for the numerical estimation of the 
central charge, in contradiction with the common belief. 



Exponential fields. We introduce a non locality factor e a i a 2 in front of the product of 
the ^-matrices in the residue equation. We propose the p-function: 

n l «2 1 1 a l "2 1 

p 2n (9,u,v) =p^ G (9,u) ®Psh 2 { e i v )+Psh( e i u )®Psh 2 ( d > v ) ( 17 ) 
+pf G +l2 (9,u) ®pf G (9,v) + P f G ~ho^)®pf G (9,v). 
One notices that the following sum of two particle form factors: 

F~ -{9 U 9 2 ) +F-T'-t(6i,e 2 ) +F~- -{9 U 9 2 ) + F^-^(8 1 ,0 2 ) 
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Pi 


c (2) 


Cexact 


P2 = 1 


1 


2 


2 




2 


1.9848 


3 


P2 = 2 


2 
3 

5.2 
15 
150 


2.1586 
2.196 
2.213 
2.217 
2.217 


3 


1)0=3 


3 
10 


2.255 
2.29 


3 


P2 = 5.2 


5.3 


2.3 


3 


P2 = 10 


10 


2.35 


3 


P2 = 15 


15 


2.36 


3 



Table 2: SS model 



has for p- function 11 

I _i I _l 
(PSG + Psh) ® (P|o + ^5g) + (PSG + Psg) ® (PSG + Psg)- 

Then, using the expression we find 

F e (0 1 ,9 2 )= (18) 

2iir 2 Nf / 7rpi 7rp 2 \ , #12 , , a \ ( (si <8> s 2 + si ® s 2 ) &> (si <8> S2 + «i <8> s 2 ) \ 
— — - — r-r cotan h cotan cosh tss["i2) 1 — ; ; 1 — ; ; — • 

{c pl c P2 y\ 2 2) 2 JssKlv \ sinh _i_ (i7r _ 0i2)sinh _i_ ( ^_ 0i2) J 

This exactly reproduces the expression for the trace operator eq. (j!6j) above, with the new 
normalization constant 



M 2 (C pi C P2 ) 4 sin^sin^ = _ q (C Pl C P2 ) 4 
(gi+gg) 77-2 



The value for the vacuum energy 



7TPl ■ 7TP2 

, sin -rr sm 



sm7r (Pl±P2) ' 



< G >= -ttM' . ' f , (19) 

2 



was first obtained by Bethe ansatz in [4]. 

From the computation made above, one may conclude that the exponential fields O ai a2 = 
e j(2ai</Ji+2a 2 </32±a3</>3) are associated to the p-function (|18|). as the expression of the trace operator 
is [4]: 

9 ^ e »«3V3 cos ( ai(/9l 4. a 2 <p 2 ) + e~ ia3<P3 cos(ai93i - a 2 ip 2 ). 



11 We used the fact (in particular here for two particles) that the form factors in the SG model with the p- function 
Psg -which corresponds to the identity operator- are equal to 0. 
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Let us consider the case p\ = P2 = 1 and P3 = 0, where we have two SG models at the 
free fermion point. Using the results of [22], the expression for form factors associated to this 
p-function gives the expected result: 



F'l 



.12 (a a \ _ COS _2£ 1 



e (v / 2(a 1 +a 2 )+i)e 21 gin ^( ai + fl2 ) _ e (V2(a 1+ a 2 )-i)e 21 sin + ^ 

= sin V2tt (at + a 2 ) e ^i+a 2 )e 21 . 

cos^f 

For more than two particles, there should be some significant simplifications in our general 
formula. 

Fateev and Lashkevich constructed in [18] with a very different method the form factors for the 
most general exponential fields of the type e 2m i¥'i+ 2ia 2^ 2 +«(2a3±a3)¥'3 ) w hich supposes that we 
introduce in the p- function a dependence with respect to an extra parameter a^, that does not 
affect the non locality factor. For the time being, I do not see how to do this for 03 arbitrary. 
The first proposal that one could imagine for the exponential fields e 2i ( ai 'P 1 + a2 V , 2±2 a 3<P3) ; w here 
k is an odd integer, is the following p-function: 

p 2n {0,U,v) = P g 1 G (e,u) ®Psg 2 (6i v )+Psh( e i u )®PsG * (.°' v ) 

+pf G +Li (0, u) ® pf G (0, v) + pf G ~ 2 (6, u) ® pf G (6, v) (20) 

However if one considers the free fermion point, one does not recover the expected result (unless 
k = 1), so this proposal should be rejected. 

It is unclear to me what the form factors of exponential fields are for 03 7^ ±203. The results 
of [18] indicate it should be possible to construct them, but I am not capable of making any 
comparisons with the results of this paper: in particular the authors of [18] have not succeeded 
in their article into making explicit evaluations of their two-particles integrals for the simpler 
cases ai = yctj, so even their two particle form factor for the trace operator is not explicitly 
computed. 

One may think of introducing a dependence w.r.t 03 by considering linear combination of func- 
tions of the type (|15|) where we would change the number of elements of T. However, such a 
function (|15|) is very close to a similar function obtained in closely related mathematical problems 
though different physical contexts in [30,31]; in these articles only the function with #T = n — 1 
was considered, and it seems to me that this is what should be done also here 12 . 



RSOS restrictions. 

• For integer values of pi and P2 > 3, the QFT admits a quantum group restriction with 
respect to the symmetry group U q2 (sl2), Q2 = exp( 2 ^ L ), giving the parafermionic sine- 

12 It is not even clear to the author whether the form factors constructed above really correspond to the 
exponential fields of the UV CFT e 2iai vi+-iia 2lp2 +±i a3V3 although we have little doubt about the trace operator- 
and are not only a formal mathematical solution of the equations. In other words, the author has no clear 
understanding of the correspondence between space of states in the CFT and space of states in the corresponding 
integrable perturbed model, beside the sine-Gordon case. 
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Gordon models [19] with action: 



(21) 



where A p 2 is the action of the Z P2 _ 2 parafermionic CFT with central charge c = 2 — 6/p 2 , 
and the fields ip, ip* , (rp,^) are the holomorphic (antiholomorphic) parafermionic currents 
with spin A = 1 — p ^_ 2 (A = —A). The field (p is a scalar boson field and the parameter 
p is given by 

2 = 8tt pi 

P ( p2 -2)(p l +p 2 -2Y 

We denote the QFT (|21|) by V(pi,p2)- The perturbing operator has conformal dimension 
Ap er t = |^ + 1 — p2 ^ 2 , In addition to the conformal symmetry, the Z p2 _ 2 parafermionic 
CFT possesses an additional symmetry generated by the parafermionic currents ip, (ip). 
Basic fields in this CFT are the order parameters aj, j = 0, 1, . . . ,p 2 — 3 with conformal 
dimension: 

s _ j(P2~ j - 2) 
3 2p 2 (p 2 - 2) ' 

All other fields in this CFT can be obtained from the fields aj by application of the gen- 
erators of the parafermionic symmetry [34]. We introduce as basic operators in V(pi,P2) 
the local fields: 

$ ( a j) = o-jexp(iatp). (22) 

The vacuum structure, spectrum and scattering theory of the QFT 1)21(1 are described 
in [19]. Its S'-matrix is 

-s^(e 12 )®s* sos (e 12 ), Pl >i. 

For p 2 = 3, ip = 1 and the QFT V(pi,3) is the usual Sine-Gordon model with coupling 
constant p 2 = 1 = §j^. 

For p 2 = 4, the parafermionic current i^{z) is a Majorana fermion and the theory "P(pi,4) 
describes the N=l supersymmetric sine-Gordon model with a coupling constant p 2 = 

Sir pi 



2(2+pi)' 

In order to obtain form factors of the trace of the energy momentum tensor in this family 
of models directly from the formula for form factors of the trace in the S'S'-model, we 
use the RSOS procedure inspired from the SG case, and we propose as p-function of the 
operator #"eW + ^t^t e -w : 

i 

pf n (6,u,v) = p%{9,u) ®Ps G (8,v), 

with the modification of the Bethe ansatz state ^efea.-.ean like i n equation ( |7| ). The two- 
particles form factor for the trace operator reads explicitly: 

^(^H^cosh^M 

P1P2 2 



(si <g> s 2 + si <8> s 2 ) <8> (e 2 "2 si (g> s 2 + e 2p 2 si (8) s 2 ) 
sinh ^-(ivr - 6>i 2 ) sinh - 12 ) 



(23) 
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It is normalized such that 

F e (9 1 + in, Ox) = 2vrM 2 (si <g> s 2 + si ® s 2 ) ® (e^ si0s 2 + e"^ si (8) s 2 ). 

The expression (|23|) for the two particles form factor of the trace of the energy momen- 
tum tensor correctly reproduces in a more or less accurate way the known value for the 
UV central charge c = 3 — 6/p2- numerical checks are presented in Table 3. One sees 
that increasing p2 with p\ fixed, the conformal dimension of the perturbing operator gets 
closer and closer to one, so one may expect a badly convergent integral ©, which would 
be responsible for the observed lack of precision in the numerical tests. For p 2 fixed and 
increasing p\, the conformal dimension again gets close to one: we do not have any expla- 
nation for the fact that for p 2 ^ 3, the accuracy becomes better for p\ bigger. 



Again, if we use the expression for PgQ, we find the following value for the normalization 



constant N®: 



~ e = M 2 (C pi C P2 ) 4 tan ^ 



np\P2 2 

It is surprising to note that had we identified the normalization of the trace operator as 
we previously did: 

n® = - < e > %^l>! 

TT Z PIP2 

this would lead to the following prediction for < >: 

<e>=-vrM 2 tan^|i, 

which is certainly true for p2 = 3, 5, 7, . . . , but not for p 2 = 4, 6, 8, . . . , where in the latter 
case the v.e.v of < O > should be vanishing, according to formula ()19|) ! Consequently, in 
this case, the normalization of the two particles form factor of the trace operator does not 
reproduce the correct answer for the vacuum expectation value of the trace operator. 

Although it might seem a priori most natural to propose for p-function for the operators 
ip'ip exp ia(p(x) 13 (and exp ia(p(x)) : 

a l 

P2n(0, U, V) = Pg G (0, U) (8) Pg G (0, v), 

with the modified Bethe ansatz state ^eie 2 ...e 2n as defined as in equation ((7J), we are not 
quite certain that this identification is correct: the problem not being a technical one, but 
a problem of existence of such operators 14 . 



We certainly expect this to be true for p = 3, in which case tp,ip = 1, and we should recover the form factors 
for the exponential fields in the sine-Gordon model. 

14 The author has already experienced a similar problem within the context of the massless N—l super sinh- 
Gordon model [33], where she came to the conclusion that it might be meaningless to consider form factors of 
e°" p ,ae°" p for a generic. 
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7^ 

L num 



Pi 



Ccxact 



P2 



1 
2 

3.6 

4 

5.2 
20 



1 

0.9924 
0.9881 
0.9878 
0.9873 
0.9869 



1 



P2 



1 

2 

3.6 

4 

5.2 
20 



1.3976 
1.4404 
1.4469 
1.4473 
1.4477 
1.4480 



1.5 



P2 = 5 



1 

2 

3.6 

4 

5.2 
20 



1.5945 
1.6819 
1.7016 
1.7029 
1.7048 
1.706 



1.8 



p 2 = 20 



1 

2 

3.6 

4 

5.2 
20 



1.9488 
2.17 
2.26 
2.27 
2.29 
2.3 



2.7 



Table 3: Parafermionic Sine-Gordon model 



P2 


c (2) 


Qjxact 


3 


0.9924 


1 


3.5 


1.2542 


1.2857. . . 


4 


1.4404 


1.5 


7.3 


1.9422 


2.1780... 


10 


2.063 


2.4 


50.3 


2.211 


2.8807. . . 


+oo 


2.217 


3 



Table 4: N=2 supersymmetric restricted Sine-Gordon model 
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Let us mention that another interesting case is obtained for p\ = 2 and P2 arbitrary integer. 
The QFT V(2,p2) possesses N=2 supersymmetry and is known as the N=2 super symmetric 
restricted sine-Gordon model. The UV central charge is the same as before. The numerical 
checks on the central charge can be found in Table 4. We have found interesting to present 
in the table above three cases where P2 is not an integer: as one may expect, the observables 
depend continuously on the parameter p 2 - 

For integer values of p\, the previous QFT T , (pi,p2) admits an additional restriction with 
respect to the second quantum group 

U qi (sh), Qi = exp(ff). The resulting QFT [19], 
which corresponds to integrable perturbations with an operator <I>c of conformal dimension 
^ = 1 ~~ P1+P2-2 °f * ne cose t m odel [20] su(2) pi _2 <8> su(2) P2 _2/su(2)p 1+P2 _4, is denoted 
by <£{pi,p2)- The S-matrix is [19,35] 



12 J 



with pi,P2 integers, pi,^2 > 3. 

For pi = 3, the QFT £(3,j>2) describes minimal CFT models perturbed by the operator 
$1,3 [25]. 

For pi = 4, the QFT <t(4,p2) describes N=l supersymmetric minimal CFT models per- 
turbed by the supersymmetry preserving operator &c- hi particular, the model £(4,4) 
possesses N=2 supersymmetry. The corresponding CFT has central charge c = 1. It de- 
scribes the special point at the critical line of the Ashkin- Teller (Z4) model. At the N=2 
supersymmetric point the thermal operator e coincides with the operator $c. 
The basic fields in the coset CFT can be represented in terms of fields aj from the 1* P2 -2 
CFT and a scalar bosonic field ip as [36]: 

$g) = a j e iatkV , (24) 

where: 

(1 — I) 9 8-7T pi 



(p 2 -2)p' r ( p2 -2)(p 1+P 2-2)- 
and the integers k < pi and I < p\ + P2 — 2 satisfy the condition 

\l - k - n{p 2 - 2)| = j; n £ Z P2 _ 2 - 

We propose as p-function of the trace of the stress energy tensor: 

/ 2n n \ , 



pf n (6, u, v) = - [J2 e~* E eUi ® PsgV> v )> C 25 ) 
\t=i i=i / 

then both Bethe ansatz states ^P 1 ^ 2 should be modified as in equation (J7J). Explicitly for 
two particles: 

F (^ 2 ) = ^cosh^/ ss( M 
P1P2 2 



x 



(e 2 Pi si (g) s 2 + e 2 w si <g> s 2 ) <g> (e 2 ? 2 s\ ® s 2 + e 2p2 s a ® s 2 ) 
sinh^-(z7r - 6'i 2 )sinh^(i7r - 6»i 2 ) 
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IT 6 






Pi 


= 3 


3 


0.5 


0.5 






4 


fiQ88 


7 






o 


D 7Q79 
u. i y i Zi 


u.o 






20 


0.9744 


0.9857. 


Pi 


= 4 


4 


0.9924 


1 






5 


1.1429 


1.1571. . . 






20 


1.4268 


1.4727. . . 


Pi 


= 5 


5 


1.3240 


1.35 


Pi 


= 20 


20 


2.2 


2.5578. . . 



Table 5: Coset model su(2) Pl _2 <8> su(2)p 2 -2/ s«(2) Pl+P2 _4 



It is normalized such that 

F H (6*i + in,0i) = 2-nM 2 (e 2p i si <g> s 2 + e 2p i si 8) s 2 ) <8> (e 2p 2 S \ <g> s 2 + e 2p 2 si®s 2 ). 
Let us note that we certainly expect for symmetry reasons the p-function 

1 I 2n n \ 

p® (9, u, v) = -p| G (0, u) ® J] e ~ 6i E eV ' ( 26 ) 

\i=i i=i / 

to lead to the same form factors as the p-function in eq. (|25jl. Also we would expect the 
p-functions 

P f n (e,u,v) = pl G {e,u)® P l SG {e,v), (27) 
P f n (e,u,v) = p 1 SG (e,u)0 P l G (e,v), (28) 

to give the same result for the form factors. 

The numerical checks provide results for the central charge in Table 5, to be compared 
with the exact result c = 3 — 6(p^ 1 +p 2 1 ~~ (Pi +P2 — 2) _1 ). For p% = 3, they are identical to 
those of Table 1. Again, increasing pi,p 2 make the conformal dimension of the perturbing 
operator closer and closer to one, leading to a badly converging integral Q, which may be 
responsible for the decrease of accuracy of the numerical results. Let us note again that 
the normalization constants iV® do not reproduce the correct v.e.v for the trace: indeed, 
with the p- function (|27|) , we would obtain 

< 9 >= -vrM 2 tan^p 



and with the p-function ()28|) 



<e>=-vrM 2 tan^, 



whereas the correct value for < > is given by (|19[) . 
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k + 2 






Q 
O 




1 
1 


A 
^± 




1 ^ 


5 


1.706 


1.8 


10 


2.1 


2.4 


100 


2.3 


2.94 


+00 


2.3 


3 



Table 6: SU(2) k WZNW model 

Certainly, if p 2 = 3 we recover the form-factors of the operators e taikip<yX ^ in the minimal 
model M pi : 

a ik i 

p 2n (0,U,v) =P S G (6,u)®p 2 SG (0,v), 

with the modification of the two Bethe ansatz states iJ/P 1 ^ 2 . 

• A common limiting case of these two RSOS restrictions is the Polyakov-Wiegmann model 
[8], which is the same as the perturbation of the SU{2) WZNW model at level k with 
action 

S = Swznw.- + A / d 2 xJ a J a , 



where J a (J a ) are the holomorphic (antiholomorphic) SU{2) currents. The S-matrix of 
the model is given by: 

-s s + g(6 12 )®s™° s (e 12 ). 



We get the Principal Chiral Field model [8,9] in the extreme limit k = +oo. Whether one 
considers S+^ s '(#12) or S+^,(#i2) amounts to the same result. The numerical estimations 
for the central charge with two particles contribution are given in Table 6. 
The Polyakov-Wiegmann model is asymptotically free, and the trace operator is of con- 
formal dimension one for any k; from the results above we deduce that the approximation 
with two particles leads to a conformal dimension of the trace operator which is less than 
one in the case of the SU(2) Thirring model (k = 1), whereas when one increases k, it 
becomes closer and closer to one already at the level of two particles. Numerically, the 
value for the central charge is given by an integral that becomes logarithmically divergent 
when the conformal dimension of the perturbing operator is one. It is not clear why for 
small k the results are quite good, and poor for big k. 



4 Conclusion 

A few remarks should be made concerning the formula (|15|) . which is one important result 
obtained in this paper. 

• The function Ai is universal as it does not depend on the scattering of the theory, and ap- 
pears in models where the S'-matrix is given by a tensor product of different S'-matrix. Typ- 
ically, such a construction (with suitable modifications) is encountered in integrable mass- 
less flows [30,33,37], which contain three different types of scattering, between right /right 
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movers, left/left movers, and right/left movers. It has also been used in a different physical 
context in [31]. 

• Other mathematical solutions are obtained by changing in formula (|15|) the number of 
elements of the set T. It is not clear to what they could correspond. 

• Of course it is crucial to make a comparison with the results of [18]. 

• The so-called one to one correspondence between operators of the UV CFT and operators 
in the massive/massless theory is a issue which deserves, to the opinion of the author, 
better understanding. 
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